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Abstract 

We illustrate the basic notions of additional non-isospectral symmetries and their interplay with 
the discrete Darboux-Backlund transformations of integrable systems at the instance of constrained 
Kadomtsev-Petviashvili (cKP ) integrable hierarchies. As a main application we present the solu- 
tion of discrete multi-matrix string models in terms of Wronskian r-functions of graded SL{m, 1) 
cKP hierarchies. 

1. Introduction. KP as an Arch- Type Integrable System 

Integrable systems (for the basics, see refs. (|, ||, ^) constitute an outstanding branch of theoretical physics 
since they describe a vast variety of fundamental non-perturbative phenomena ranging from D = 2 (space-time 
dimensional) nonlinear soliton physics and planar statistical mechanics to string and membrane theories in 
high-energy elementary particle physics. It turns out that, under plausible assumptions, a variety of physically 
interesting theories in higher space-time dimensions can be reformulated as lower-dimensional {D = 2) integrable 
models which in the same time possess infinite- dimensional symmetries and thus, as a rule, being integrable 
(see, especially, the recent developments Q related with intcgrability of Seiberg-Witten effective low-energy 
theory of (extended) supersymmetric gauge theories). 

Among the various infinite-dimensional symmetry groups and algebras playing role in integrable field theo- 
ries, a particularly distinguished place belongs to the Lie algebra Wi+oo (|] (specific "large N limit" of Zamolod- 
chikov's Wn conformal algebras isomorphic to the Lie algebra of all purely differential operators on the 
circle). It contains (together with its supersymmetric extension) all previously known infinite-dimensional 
symmetry algebras - Virasoro and Kac-Moody. Also, it is precisely a Lie-algebraic deformation of the infinite- 
dimensional generalization of the Virasoro algebra - the algebra of area-preserving diffeomorphisms. 

Recently Wi^oc symmetries attracted broad interest as they appeared naturally as inherent structures 
of models in different areas of theoretical physics : theory of black holes and space-time singularities, two- 
dimensional quantum gravity, nonlinear evolution equations in higher dimensions, self-dual gravity, N = 2 
superstring theory (refs. (a)-(d) in 0), quantum Hall effect All listed models possess, in one form or in 
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another, exactly soluble features, which naturaUy suggest an intimate connection of integrability with VFi+oo 
algebra. This claim may, furthermore, be substantiated by the realization that VFi+oo algebra is a subalgebra of 
the algebra '^VO of arbitrary pseudo-differential operators on the circle. Already since the pioneering papers 
of Adler-Kostant-Symes and of the Faddeev's school |^ it was realized that "^VO forms the foundation of 
completely integrable systems. In fact, it turns out that most known integrable models (i.e. those admitting 
Lax or "zero curvature" representation) can be associated with specific coadjoint orbits of various subalgebras 
of "^VO or with different Hamiltonian reductions thereof |Q . 

The generic integrable system based on "^VO symmetry algebra is the Kadomtsev-Petviashvili (KP) in- 
tegrable hierarchy |ll| of soliton nonlinear evolution equations. Its name derives from the fact that KP 
hierarchy contains the D = 2 + 1 dimensional nonlinear soliton KP equation which appeared originally in plasma 
physics. In the last few years the main interest towards KP hierarchy originates from its deep connection with 
the statistical-mechanical models of random matrices ((multi-)matrix models) providing non-perturbative dis- 
cretized formulation of string theory [p^ . 

The purpose of the present talk is to provide a brief discussion of the basic notions of additional non- 
isospectral symmetries |l^, ^ and their interplay with the discrete Darboux-Bdcklund transformations JT4[ 
of integrable systems within the KP integrable hierarchy (being and arch-type integrable system, as pointed 
out above) and its various constrained versions (cKP hierarchies) relevant in discrete multi-matrix string mod- 
els. Furthermore, we show how to obtain the solution of the latter in terms of Wronskian r-functions of the 
cKP hierarchies. 

2. KP Hierarchy: Pseudo-Differential Operator Formalism and Hamiltonian 
Structures (W-Algebras) 

We describe the (generalized) KP integrable hierarchy in the language of pseudo-differential operators (for a 
background, see (ll)). The main object is the pseudo-differential Lax operator L subject to an infinite set of 
evolution equations: 

m— 2 „ 

j"=0 i>l ' 

Here, the coefficients of L are (smooth) functions oi x = ti and the higher time-evolution parameters ^2, ^3, ■ ■ • 
] D = dx, whereas the subscripts (±) denote purely differential (purely pseudo-differential) part of the cor- 
responding pseudo-differential operators. The flows ^ in ([!]) commute (as vector fields on the space of Lax 
operators (|l|)) among themselves which expresses the integrability of the KP system. 

Within the Sato- Wilson dressing operator formalism, with the following dressing expression for the gener- 
alized KP Lax operator (|^): 

L^WD'^W-^ , W = l + '^w.,D-' (2) 



i>i 



the evolution equations for L are equivalent to: 



= -(WD^W^^) W (3) 
dti V V ^ 

In what follows we shall also need the important notions of (adjoint) eigenfunctions and (adjoint) Baker- 
Akhiezer functions. The function <I> (^P) is called (adjoint) eigenhmction for the Lax operator L satisfying Sato's 
flow equations (|^) if its flows are given by expression: 

for the infinite many times ti ^. If, in addition, an (adjoint) eigenfunction satisfies the spectral equation: 

^-V'(A) - L^^jV^A) ; L>*(A) =AV*(A) , 
'0'^*)(A) is called (adjoint) Baker- Akhiezer (BA) function 



LV'(A) = AV'(A) , — V'(A) - jV^A) ; L>*(A) = AV*(A) , g^rW = - (i*)(+) rW (5) 



Here and below, the superscript "*" on operators indicates pseudo-differential operator conjugation. 
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The BA function of the generic Lax operator L is obtained from the BA function of the "free" Lax 
operator L(°) = : 

V^W(A)=cxp{^iA-}^e«"*^^^^ (6) 



i>i 



by applying the dressing operator W (||) to (||) 



^A) we ^^^^^^^ 



(7) 



The function T{{ti}) of aU evolution parameters is called r- function of the (generalized) KP hierarchy and by 
itself constitutes an alternative natural way to describe the pertinent integrable system. 

The KP system (^ is endowed with bi-Hamiltonian Poisson bracket structures (another expression of its 
integrability) which results from the two compatible Hamiltonian structures on the algebra of pseudo-differential 
operators ^VO 113]. The latter are given by: 



{{L\X),{L\Y)}^ = -{L\[X, Y]) 
{ {L\ X) , (L| Y) }^ = Tr 



YL 



(LX)(+) LY - (XL)(+) 
dxYies{[L, X'\jd-'^Res{[L , Y]^ 



(8) 



(9) 



where the following notations are used. < -j- > denotes the standard bilinear pairing in "^VO via the Adler 
trace {L\ X) = Tr^ {LX) with TtaX = jResX. Here X^Y are arbitrary elements of the algebra of pseudo- 
differential operators of the form X — X)fc>-oo ^''^k and similarly for Y. The second term on the r.h.s. of (||) 
is a Dirac bracket term originating from the second-class Hamiltonian constraint Vm-i = on L (|l|). 

In terms of the Lax coefficient functions Wm_2, • ■ ■ , "foi ""ij ^2, . . . , the first Poisson bracket structure (H) takes 
the form of an infinite-dimensional Lie algebra which is a direct sum of two subalgebras spanned by {vj} and 
{ui}, respectively. The latter is called Wi+oo-algebra [^. Its Cartan subalgebra contains the infinite set of 
(Poisson-)commuting KP integrals of motion — jTr^L~ whose densities are expressed in terms of the 

T-function (0) as: 

9,J-lnr = ResL^ (10) 
Oh 

In turn, the second Poisson bracket structure (||) spans a nonlinear (quadratic) algebra called Woo('ti) |]16| , 
which is an infinite-dimensional generalization of Zamolodchikov's conformal algebras |6j. 

3. Constrained KP Hierarchies. Free-Field Realizations 

Let us now turn our attention to a specific class of Hamiltonian reductions of the full (generalized) KP system (P 
(cKP hierarchies, for short), where the purely pseudo-differential part of the KP Lax operator is parametrized 
through a finite number of functions (fields). To this end let us recall the notion of (adjoint) eigenfunction of a 
Lax operator (^. As in let us consider the fiow of a vector field da given by: 



M 



daL 



(11) 



where $i, are a set of AI independent (adjoint) eigenfunctions of i (||). Using the simple identity valid for 
any differential operator i?(+): 



one can easily show that: 



9q, 



dt. 







/ = 1,2, 



(12) 



(13) 



Now, the constrained KP hierarchy (denoted as cKP^^m) is obtained by identifying the "ghost" fiow da with 
the isospectral fiow which, upon comparison of (|l^) with (|l|), implies the following constrained form of L: 



M 



L = L. 



m,M 



= L 



(+) 



m-2 



3=0 i=l 



(14) 
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subject to the same Lax evolution equations as in (^. Moreover, using again identity (12) one finds that the 
functions <i>i, remain (adjoint) cigcnfunctions of the constrained Lax operator Lm,M (14). 

As shown in ref.llj], the cKPm,M hierarchies given by dl4) are equivalent to the so called "multi-boson" 



cKP hierarchies |19 



L 



(+) 



M 



D-B 



{M)\ 1 



M, 



where 



B 



, fk] = det 



s=l 

—dx In 



M, 



fc+lj 



(15) 

(16) 
(17) 

(18) 



denotes the standard Wronskian. 

There is still another useful representation of the cKPm,M Lax operator as a ratio of two purely differential 
Lax operators [Eol |ll|2|: 



.,AI 



L 



m+M 



m,M >1 



(19) 

(20) 



L„,+M ^{D- h„,+M) {D - b,n+M-l) ■■■{D-bi) , Lm = [D- bM ][D- bM-l ]---[D-bi\ 
where the coefficients bj^bj are subject to the constraint: 

m+M 

Y^b,-Y,h^Q (21) 

3 = 1 1=1 

As already proved in detail in re f . pl| ] , the cKP^.m Lax operator L — L.^.m obeys the same two compatible 
Poisson bracket structures ^ and (|9|), i.e., cKP„i j\/ hierarchies are legitimate Hamiltonian reductions of the 
full (generalized) KP hierarchy. Moreover, the second Poisson bracket structure (^ in terms of the coefhcients 
(pO|) takes the form of free-field Poisson bracket algebra: 



{ h{x) , bj{y) } = - ^ 6'{x - y) , i,j 
{bu{x),bi{y)] = -[5ki + ^ S'{x~y), 
{h(x),bi{y)] ^ -5\x-y) 



= l,...,m + M 
kJ^l,...,M 



(22) 



which, as demonstrated in refs.|g3|, |2l|, is precisely the Cartan subalgebra of the graded SL{m -f M, M) Kac- 
Moody algebra. This latter property justifies the alternative name of the constrained cKP„jjvf hierarchies - 
SL[m + M, M) KP-KdV hierarchies. 

In other words, (p^-(p2|) provide via eq.(^ explicit free-field realizations of the nonlinear 'W^{m) algebra. 
Similar free-field realizations exist also for Wi+oo ^ the first KP Poisson bracket structure (see refs.[^ pT|). 

4. Additional Symmetries and Darboux-Backlund Transformations 

Let L be again a pseudo-differential Lax operator of the full generalized KP hierarchy (|l|) (recall x = ti) and 
let M be a pseudo-differential operator "canonically conjugated" to L such that: 



L, M 



dti 



^(+) ' M 



(23) 



Within the Sato- Wilson dressing operator formalism (||)-(||) the M-operator can be expressed in terms of 
dressing of the "bare" M'"' operator: 



z — / rfji ^ — / 
i>i i>i 



X. 



I 



(m) 



(24) 
(25) 



1=1 
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conjugated to the "bare" Lax operator L'^*'^ — D™- , i.e.: 



m 



m I \ ^ I + m - 



l>0 



m 



m dti 



(26) 
(27) 



where m ( |27| ) we used eqs.(|3|). Note that X(,„) is a pseudo-differential operator satisfying 
On BA functions (lEI) the action of M is as follows: 



(m) 



= 1 



MV'(A)= (^|^ + «™(A)) V(A) (28) 

where am{\) is a function of A only^. 

Since any eigenfunction <i> (^) can be represented as a linear "superposition" of BA functions (^) : 

$({n)=^dA^(A)V'(A,{t}) (29) 



The so called additional (non-isospectral) symmetries |ll| are defined as vector fields on the space of 
KP Lax operators (0) or, alternatively, on the dressing operator (0), through their flows as follows: 



(with an appropriate contour F in the complex A-plane, such that the integral in ( p9[ ) exists), eq.(|28|) implies 
that: 

M^{{t}) ^ j^dX (^-^ + a„(A)^ 0(A)VXA, {t}) (30) 



dk.nL = — 



'(+) 



k 71 ,f n — 1 



nL^M 



(31) 



which commute with the usual KP flows ^ (|^) . 

Let us now turn our attention to the notion of Darboux-Backlund (DB) transformations of (generalized) 
KP hierarchy ([|) and its reductions - cKP„i_j\/ hierarchies defined as follows jlj, |l^ : 



L 



TLT^^ = L 



(+) 



L_ 



T = xDx' 
W+x{d,{x-'h+)X)^,D-^) X-' 



(32) 
(33) 



<i>oD-'^o+xDx-'L^xD-'x-^ = <i>oD-'^o + J2^,D-^^', for i = L,„,m j 



$0 = (x^X^'i) X = TLX 

^^ = Xa.(x"'l'0 



*o = X ^ 
-X-'d-' (*,X) 



(34) 

(35) 
(36) 



where x is an (non-BA) eigenfunction of L. The DB-transformed Lax operator (p3) satisfies the same flow 



9ti 



equations w.r.t. ti as in (|1| 

operator i? 

(xi?x"'Sx^"'x"' 



due to the simple identity valid for any pseudo-differential 



(+) 



XDX-'B^+^XD-^X-' - Xd. (x~HS(+)X)) D^^X'^ 



(37) 



(and using the fact that x is an eigenfunction of L). Moreover, eq.(p4[) shows that, in order to preserve the 
cKPm,M form of the DB-transformed Lax operator L — Lm.M, we have to choose x = where is any one 
of the eigenfunctions of the initial L = Lm.M (p^). 

One can generalize (32)-(|3^) for successive Darboux-Backlund transformations on the initial L — Lm,M = 
i*^"^ as follows. Within each subset of m successive steps we can perform the DB transformations w.r.t. the m 



^The appearance of am (A) can be traced back to the ambiguity in the definition of the dressing operator (|2|): W 
WWo where Wq = 1-1- X/i>i ^if"' with constant coefficients Ci. 
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different eigenf unctions of (^4|). Repeated use of the following important composition formula for Wronskians 
PI : 

Wk{f) 



TkTu-i ■•■ Ti(/) 



where 



£»+ In 



Wk = Wk[4'u.--,4'k] , Wk-i{f) = Wk[^Pi,...,i,k-iJ] 

and employing short-hand notations: 



is) _ 



(0) 



Z = 1, . . . , 771 



(38) 

(39) 
(40) 

(41) 



where the upper indices in brackets indicate the order of the corresponding DB step, yields the following 
generalization of (^5|)-(^6|) (below 1 < / < tti) : 



I (fcm+i) _ rp{k7n+l-l) rp(km) rp(km-l) rp((k-l)m) rp(m-l) rp(n) {k±) 
~ . ■ ■ ±1 J„ . . . •••Jm ■ ■ ■ Xi 

VY ^1 J ■ • ■ , , Xl J ■ • ■ I Xm , ■ • ■ 7 Al : ■ • ■ j Xm , Xl : ■ ■ ■ i Xl : Xi 



w 



^(0) ^(0) (1) (1) (fe-i) (fe-i) (k) (k) 

, • ■ ■ , 1 Xl J • ■ ■ I Xm , ■ • • , Al I • • ■ I Xm t Xl i ■ • ■ i AZ 



(42) 



(fc+) (fc+l) r 1 ^ • ^ ; 

Xi = Xi for 1 < 7 < / 



xl'^ ' = Xi*^^ for / + 1 < 7 < 771 



Correspondingly, for the r function (0) after km + Z steps of successive DB transformations we get: 

^(km+l) 



t(o) 
W 



,(fem+/-l) <j,(fc™)<j,(fcm-l) ^{{k-l)m) ^(m-l) <j,(0) 



^ 1 1 ■ ■ ■ 1 ' Al 1 ■ ■ ■ 7 Am 7 ■ ■ ■ 7 Al i ■ ■ ■ 7 Am ■> Al 



7 ■ ■ ■ 7 



Xi 



(fc) 



(43) 



We now formulate the main result of this section - the condition for compatibility between additional- 
symmetry flows (l3l| ) and Darboux-Backlund transformations (^). 

Let $ be an eigenfunction of L defining a Darboux-Backlund transformation, i.e. : 



d -L ~ 

dti (+) ' V ; V ; 

or, in terms of dressing operator: 

W = ($i:>$"i) W D-^ 
Then the DB-transformed M operator (cf. (p6|)) acquires the form: 

iVf = ($i:>$"i) M {<i>D-^<i>-^) = ^-^^t„,.+iLfl^ + M- 



M- - WX(,^)W-^ -im-Y] -tm+i . W-^ 

^ ' ^-^ m oti 

i>i 



(44) 
(45) 

(46) 
(47) 



where ^(m) — -DX(„i)-D^^ with X^^-) as in (p5|). Clearly Xj-^-j, like X(„i), is also admissible as canonically 
conjugated to D™. 

The DB-transformed BA function reads in accordance with (Eq): 



7/-(A) = A-^$9,($-V(A)) 
and the DB-transformed Af-operator acts on it as: 

MV^(A) = + "m(A)) ^(A) , 5„,(A) = «„,(A) + i-A"! 



(48) 



(49) 
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Taking into account (^), we arrive at the following important 
Proposition. Additional symmetry flows (^Jj commute with Darboux-Bdcklund transformations (^^-(^^, 



I.e. 



dk,nL 



dk,nW 



if and only if the DB-generating eigenfunction <f> transforms under the additional symmetries as: 



'(+) 



(50) 



(51) 



Motivated by applications to ( mult i-) matrix models (see next sections and ref . pOj ) , one can require invari- 
ance under some of the additional-symmetry flows, e.g., under the lowest one 9o,i known as "string-equation" 
constraint in the context of (multi-)matrix models: 



do,iL = 



M(+) , L 



= -1 



5oa$ = 



M(+)$ = 



(52) 



Eqs.(p2D, using second eq.(|23|),(|26D-(p9D and ( |5l[ ), lead to the following constraints for L, the BA function tp{X) 
and the DB-generating eigenfunction $ of L, respectively: 



El + m d 
m ati 



ti,L 



-1 



■sr-^ I + m d , ,\ , , d , , 

j>i ' 



\ ^ I + m d I „ 

> im+i 77- + <m I = 

1>1 ' 



$( 



({i})= / dAe/"(^V(A;W) 



(53) 
(54) 
(55) 



Now let us recall the formula ( ^ ) for the r-function ratio for cKP„i,a/ hierarchies subject to successive 
DB transformations. Noticing that the eigenfmictions <i>*^'^-' of the DB-transformcd Lax operators L*^*^^ satisfy 
the same constraint eq.(p5[) irrespective of the DB-step k, we arrive at the following result ("string-equation" 
constraint on the r-functions) : 



Proposition. The Wronskian r-functions ^4^ ) of the cKP„i.m hierarchies satisfy the constraint equat 



ion: 



\ - I + m d 



(56) 



Example: It is well-known that the discrete one-matrix model can be associated to the following chain of 
the Lax operators connected via DB transformations: 



^ $(fc) (^\j,ci>i'')y + (^c^^y ^ ^(fc+1) ^ (^^ik)y^ 



where 



(57) 
(58) 

Ty„[0,9(/.,...,9"-V] ' ^(0) - ""LV',^^^,---,^ ^> (59) 
with (j) = J dA exp (^^-^ ifeA*"') . The above proposition (with to = 1) coincides perfectly with the "string- 
equation" £L^%7v[</',50, . . . ,9^-^] = 0, with = J2T=2 '^^featflT + ^^1- 

5. Multi-Matrix Models as cKPm,i Hierarchies 



_ w^+i[ci,,d<j>,...,d"4,] 



The partition function of the multi-matrix (g-matrix) string model reads: 

9 Pa 



/( q Pa <?-i 
dMi ■ ■ ■ dMg exp - <^ ^ ^ Tr M^" + ^ Tr Af„M„+i 



(60) 
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where Ma are Hermitian N x N matrices, and the orders of the matrix "potentials" Pa may be finite or infinite. 
In refs.|Q it was shown^ that, via the method of generahzed orthogonal polynomials one associates to 
( |60| ) generalized Toda-like lattice systems subject to specific constraints, so that Zn and its derivatives w.r.t. 
the coupling parameters can be expressed in terms of solutions of the underlying Toda-like discrete integrable 
hierarchy where {i^^^}, . . . , {i^'''} play the role of "evolution" parameters. This Toda-like discrete integrable 
hierarchy differs from the full generalized Toda lattice hierarchy in that the associated Toda matrices in 
the first hierarchy are semiinfinite and contain in general finite number of non-zero diagonals. 

It turns out that, in order to identify the continuum cKP integrable hierarchy which provides the exact 
solution for (|60|), we need the following subset of the associated linear system and the corresponding Lax 
( "zero-curvature" ) representation from the Toda-like lattice system : 



d 



Otr '- 



d 



dt 



(9) 



d 



d 

—)Q{<i) 



Q\+YQ{q) 



dti 
d 



dt 



In what follows it is convenient to introduce the short-hand notations: 
^=41) , r = l,...,pi ; is =4'^ , s = l,. 



^Q(l)= [Q(l),Q(g)l] 
[Q(9),Q(<?)1] 

g = Q(l) , Q = Q{q) 



(61) 
(62) 
(63) 

(64) 



Further, there is a series of additional constraints ("coupling conditions") relating Q = Q{1) and Q = Q{q). In 
the two-matrix model case {q = 2) their explicit form is: 



Qw = - E 



s=l 

pl-1 



9 

(r + 1) 



9 

s+iQ 



Q:Q 



1 



(-) 



-ill 
9 



tr+lQ^+) ill 

9 ^ ' 9 



(65) 
(66) 

(67) 



Here the subscripts — /+ denote lower/upper triangular parts, whereas (-!-)/(—) denote upper/lower triangular 
plus diagonal parts. In the higher (g > 3) multi-matrix case the "coupling conditions" have much more intricate 
form (involving also the "intermediate" Q{2), . . . , Q{q — 1) matrices). However, their explicit form will not be 
needed to find the solution for Z^f (pO[) since we will be able to extract the relevant information solely from 



the discrete Lax system (|62 



-(|63|) and the relations expressing Q = (5(1), Q = Q{q) in terms of orthogonal 



polynomial factors (see eqs.(7|) below). 

The parametrization for the matrix elements of the Jacobi matrices Q = Q{1) and Q = Q{q) is as follows: 



Q(l)n™ = 



= Q(l)„„ = Mn) 
]n,n^k ^ a-kin) k 



QWn,n+l = Qn,n+1 - 1 

l,...,m(l) , rn(l) = (pg 
for m ~ n > 2 , n — m > m(l) + 1 

, = 6o(n) , Qiq)n n-l = Qn,n-1 = Rn 



Q(.l)n^n+k = Qn,n+k = 6fc(n)i?„+i ■ • • i?„+fe fc = 1, . . . , m((j) , m{q) = [pq^i 

Q('i)nm = Qnm = for n-m>2, m~n> m{q) + 1 



1) • . • (P2 - 1) 



1)...(P1-1) 



(68) 



(69) 



In terms of the Q — Q(l), Q — Q{q) matrix elements the partition function (pQ) is expressed in the following 



way 



N-l 



Zn = const hn 



(70) 



n=0 



"See also refs.[^ and the lecture of Prof. L.Bonora in the present volume. 
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ao{n) 



dti 



■ In hr, 



Oti hn-l 



Otr 



d_ 



N-1 



}nn I or ^" ~ Qr, 



N-1 



(71) 



— InZw = ^ ao{r 



n=0 
iV-1 



T^ln^AT = V ba{n) 
oti ^ 

n— 



9 



JV-l 



9t 



In Zn =^ Qnn . ^ = ^ Q 



n=0 



(72) 



n=0 



where /i„ are the normaUzation factors in the nonlocaUy generalized orthogonal polynomial formalism [p8[ (using 
notations (|64|)) : 



/ / dAd/.P„(A)exp{^A%}p(A,/x;{i"},{g}){^^%}p„(^) (73) 

/■J— n-1 Pa g-2 
ffa.a+li^ai^a+l + .912^(^2 + 9q-l,ql^q-lfJ- \ (74) 
- —1 - —I a=2 

{t"}= {t<^^\...,t<^'^~ 



a=2 ra=l 



As in the case of two-matrix model |2l|, |18|, using the lattice equations of motion (eqs.(62)-(|63|) for r = 
1, s = 1) we obtain the following important: 

Proposition. The matrix elements of Q = Q(l) are completely expressed in terms of the matrix elements of 
Q = Q{q) through the relations: 



''(-) 



(-) 



^f-) ' = E ^--QU ' 5 = 0, 1, ... , m(l) 



(75) 



where the coefficients q;o,7s,o are ti -independent, whereas the coefficients as,jscr with cr > 1 are independent of 
ti,ti. All 7s(j are expressed through as = 7m(i),s solely: 



Iss = iniY , 7s.s-i = s(7ii)* Sio , 7s,s-2 = (711) 



NS-2 



7ll = (Q!m(l)J 



g(g - 1) f n2 . 

(710) +s 



n(l) 



710 



731 2 
"m(l)-l 



(76) 



'^(1) (am(i)) 



731 2 

7io 



37 



to(1) 



am(i)-2 m{l) - 1 Q^m(l)-l 



K(i)) 



2m(l) 



im(l). 



™(1) 



(77) 



In the two-matrix model the explicit form of the coefficients as reads: 0:5 = —^^^^^is+i- 

Similarly, we have the dual statement with the roles of Q = Q{1) and Q = Q{q) interchanged. 

As an important consequence of ([75|), let us take its diagonal 00-part and use the last eq.(71) which yields: 

d /"^^^ 9" \ 

Of^f^o = lYl^s-^ + ao \ ho (78) 

This equation is the only remnant of the constraints ("coupling conditions") on the multi-matrix model Q- 
matrices which will be used in the sequel. 

Based on our experience with the two- matrix model ||^, Q , it turns out natural to introduce the fractional 
power of Q = Q(l) : 

g = q4t7 = Q(1)4t) (79) 
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whose parametrization closely resembles that oi Q = Q{q) m_ 



Qnn — bo{n) , Qn,n-1 — Rn , Qn,n+k — &fe("-)^„+l • ■ ■ Rji+k ^ > 1 

Qnm = for n — m > 2 



From eqs.(75) we find the following relation between the matrix elements of Q = Q(l)'"(i) and Q : 



Rn=liiRn , feo(") = 711^0 ("-) +7io , bi{n) ^ "flibi{n) + 



731 



37: 
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7io 



(80) 



(81) 



etc., with 7-coefficients as in (|7^)-(|77|). 

In order to identify the continuum cKP hierarchy associated with the general g-matrix model, as a first 
step we reexpress, using (|7|), the Toda-like lattice hierarchy (|6l|)-(|63|) as a single set of flow equations for 
Q = Q(1)4t) : 



d_ 



(QI-)) 



Q,Q- 



s = 1, . . . 2to(1), 3to(1), . . . ,pim(l) 



rm(l) 



for r = 1, . 



,Pi 



(82) 

(83) 
(84) 



Here, as in the two- matrix case |l8| , we have introduced a new subset of evolution parameters {is} instead 
of |ts = ti'^^ I defined as: 



_d_ 

dis 



d 



,m{q) 



(85) 



As a second step, one employs the Bonora-Xiong procedure |26 to get from 



with the same jsa a-s in (fi 

the discrete Lax system (B2|)-(83) an equivalent continuum Lax system associated with a fixed lattice site n, 
where the continuum space coordinate is a; = ii. Namely, the latter continuum integrable system is obtained 
by writing eqs.(p2[) in more detail using the parametrization (|6^),(|80|) : 



(86) 
(87) 



P2-1 



Xipn = Tpn+i + ao{n)ipn + ak{n)ipn-k 



fe=l 



-Rni'n 



and further using (87) to express 'ipn±e in terms of ipn at a fixed lattice site n in eq.(|86|) and the higher evolution 
eqs.(|8^) (for s > 2). Upon operator conjugation and an appropriate similarity transformation, it acquires the 
form (as before x = ti) : 



^L(n) = [(L4i)(n)) , L(n) 1 , s = 1, . . . 2m(l), 3m(l), . . . ,pim(l) 
Ots L V / (+) J 

L{n) = D^^^^ + m(l)6i(n)i?™(i)-2 ^ . . . _^ j^^^^ ^j^^ - 6o(n)) 



(88) 
(89) 



where &o,i(^),^n+i are the matrix elements of Q (p0|),(|8l|). Rewriting (|89|) in the equivalent "eigenfunction" 
form: 



L{n) = + m{l)h{n)D^^^'>-^ + • • • + $(n + l)I?^^*(n + 1) 

$(n + l) EE^„+iexp{ /6o(»^)} , + 1) = exp|- / 6o(»^)} 



(90) 
(91) 



and comparing with (|lj), we identify the continuum integrable hierarchy (|88|), describing equivalently the 
discrete multi-matrix model, as a constrained cKP ,„(!). i hierarchy. 

6. Partition Functions of Multi-Matrix Models: Darboux-Backlund Solutions | 

more detailed presentation of the material in this section will appear in ref . |3C|] . 
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Exactly as in the two-matrix case [£l|, [l^ , lattice shifts n ^ 71 + 1 in the underlying discrete Toda lattice system, 
described by (|6l|)-(pl|), generate Darboux-Backlund transformations in the continuum cKP m(i),i hierarchy 
(p8|)-(p9|). This is due to the fact that the latter continuum hierarchy preserves its form for any value of the 
discrete label n. The solutions for the eigenfunctions and r-functions at each successive step of Darboux- 
Backlund transformation is given explicitly, as particular cases of eqs.(^)-(^3|), by: 



$(n) 



Wn+i [<&(o), i^(-i)$(o), . . . , L{-ir^m 

Wn [$(0), L(-l)<i>(0), . . . , L(-l)"-i$(0)] 



|-$(o) = {L{~i))ii-^ m 



= 1, . . . ,P2,2(P2 - 1),3(P2 - 1),. . .,Pl(p2 - 1) 



(92) 
(93) 
(94) 



where everything is expressed in terms of the eigenfunction $(0) of the "initial" Lax operator L(— 1). The 
difference with the two-matrix case is only the explicit form of the latter (recall x = ti) : 



-710*1 



(95) 



where the coefficients as,7iO:7ii have more complicated dependence on {ts\ than in the two-matrix case. 

Exactly as in the two-matrix case, we obtain the relation between the n-th step DB eigenfunction $(n) and 
the orthogonal polynomial normalization factor hn (73) which generalizes (p6|) : 



$(n) = e/ = K 7n exp {£1710 + e{i')] , t' = (^2, ■ • • 
Substituting (|96|) into (|70|) and using the Wronskian formula (|93|) we get: 



N-l 



Zn=\{ =det —-^[L{-l)y~'m 



dil 



-A^(ti7io-(-e(t')) 



7ll 



det 



__ (e-*^^-L(-l)e*i^-) K 
dti ^ ' 



where we absorbed the 711-factors via changing ^ * ^ ^^^^ definition (pSf), i.e., 7^^^ ~ 'ST 

find using (Isl) and (|7|) : 



^m(l) 



J-1 



s=0 



fin — ■■ — T iT'Ci 



(96) 

(97) 
(98) 

Now, we 

(99) 



Substituting ( |99| ) into ( |98| ) yields the final result for the multi-matrix model partition function: 



det 



(100) 



which is functionally the same as for the two-matrix model, however, with a more complicated expression for 
^0 dZl) : 

Pl Pq 



o^J J dXdti exp{|^A%}p(A,/.;{t"},{g}){X]M%} 

^ ^ r— 1 s— 1 



(101) 



Eq.(lOO) was previously obtained (see refs.lSlI) from a different approach. 
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